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Introduction 


Radiation heat transfer in combustion systems has been receiving increasing interest (ref. l).In the case of hydrocarbon 
fuels, a significant portion of the radiation comes from soot particles, justifying the need for detailed soot formation 
model and radiation transfer calculations. For laminar gas jet diffusion flames, results from this project (4/1/91- 
8/22/95) and another NASA study (ref. 2) show that flame shape, soot concentration, and radiation heat fluxes are 
substantially different under microgravity conditions. 

Our emphasis is on including detailed soot transport models and a detailed solution for radiation heat transfer, and 
on coupling them with the flame structure calculations. In this paper, we will discuss the following three specific areas: 

1. Comparing two existing soot formation models, and identifying possible improvements. 

2. A simple yet reasonably accurate approach to calculating total radiative properties and/ or fluxes over the 
spectral range. 

3. Investigating the convergence of iterations between the flame structure solver and the radiation heat 
transfer solver. 


Modeling^of Tet Diffusion Flame Structure 


The structure of turbulent jet diffusion flames is modeled using the Favre-averaged equations for conservation of 
mass, momentum, and mixture fraction. A conserved scalar approach (refs. 3, 4 and 5) with an assumed probability 
density function (pdf) and a k-e-g turbulence model (ref. 6) are used. All governing equations can be written in a 
general form as (ref. 7) 
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where <f> = 1, u (velocity), f (mixture fraction), k (kinetic energy), e (dissipation), or g (variance on/). Details for 
(effective viscosity), (source term) and assumptions can be found in references. Buoyancy effects are considered in 
mean flow only, neglecting buoyancy-turbulence interactions. For laminar flames, Eq. (1) is simplified accordingly. 


The system represented by Eq. (1) is solved using the block-tridiagonal code of Chen et al. (ref. 8). State relationships 
are constructed from equilibrium calculations using STANJAN (ref. 9) to allow chemical reactions to be decoupled 
from flow calculations. More accurate laminar flamelet approaches (refs. 10 and 11) may be considered in the future. 
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We have considered two sets of soot formation and oxidation models. Both describe the transport of soot particles 
based on Eq. (1). The two-equation model developed by Moss, Syed and Stewart (refs. 12 and 13) is based on number 
density (N) and volume fraction (f v ), and the two respective source terms are 
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where = 6xl0 26 is Avogadro's number and p, is the mass density of soot (typically 1. 8-2.0 g/cm 3 ). Rate-equations of 
Arrhenius type are used to model nucleation (a), growth ( P), coagulation (y) and oxidation (x), in terms of flame 
temperature T, mixture density p , fuel and gas mole fractions, and activation temperatures. Existing soot oxidation 
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models (refs. 14 and 15) are modified to include oxidation by both 0 2 and OH. Details on model coefficients and the 
oxidation model can be found in ref. 16. For laminar flames, both Eq. (1) and rate-equation models are simplified, 
with the addition of a radial thermophoretic velocity v, =-.0S4(vlT)(oTldr). 


Another model is a one-equation model, based on volume fraction only, developed by Khan et al. (refs. 17 and 18). 
This model characterizes soot formation by an Arrhenius-type equation, and the corresponding source term is 




(4) 


where C* = 16.8 kg/Nms, Pju denotes the partial pressure of unbumed fuel, <p is the local fuel/ air equivalence ratio, 
and E s = 40000 cal/mole is the activation energy. The soot oxidation model of Lee et al. (ref. 19) is adopted, which is 

similar to Eq. (4), except with <p 3 (plp s ) replaced by / v /V? . 


The model by Khan et al. is simpler, but it may not be as physically sound as that by Moss et al., since it does not 
include such mechanisms as coagulation which causes decreasing number density under constant volume fraction. 


Solution fo r Flame Radiation Heat Transfer 

Although it plays an important role, radiation heat transfer is usually not treated in detail in combustion analyses 
due to such difficulties as the accuracy and computational efficiency of various solutions for the radiative transfer 
equation (ref. 20), the coupling of radiation and combustion solvers, the accuracy of spectral radiative properties, and 
efficient methods for integrating fluxes and other results over the spectral range. 

We chose the recently developed YIX method (refs. 21 and 22) for calculating the radiative heat flux and its 
divergence. The YIX method is a numerical approach for solving the integral formulation of the radiative transfer 
equation by reducing the order of the multiple distance-angular integrals. The name YIX comes from the shape of the 
pattern of integration points for three, two, and four angular directions. One important attribute of this method is 
that these integration points can be pre-calculated and stored, significantly reducing computational time. For multi- 
dimensional geometries, discrete ordinate sets are used for angular quadratures. Although computationally quite 
intensive, the YIX method has proven to be very accurate and suitable for nonhomogeneous media (refs. 23 and 24). 
Simplified models such as a fixed percentage of local heat loss from equilibrium everywhere, a temperature modified 

from that under equilibrium as T = T eq [l- p(T eq I T eqmax )*] (ref. 14), and a heat '‘sink" as 5^ = -eoiT 4 - Tt )/t (ref. 25) 
have each been tested and found inaccurate. We derived the spherical harmonics (Pn) approximate solution for non- 
homogeneous media following ref. 26, but found it numerically unstable in optically-thin regions. 

The radiation solver is coupled to the flame structure solver through the energy equation. For diffusion flames, the 
energy equation takes the same form as Eq. (1), with ^ = H (total enthalpy) and S ^ = —V • q rad , and the latter can be calcu- 
lated from a solution for the radiative transfer equation. These two solvers, both depending on the temperature, are 
computationally incompatible. The solver for flame structure, soot transport, and the energy equation solves 
parabolic differential equations of Eq. (1) type. The solver for radiative transfer, YIX or otherwise, is inherently not 
parabolic, since radiation is an "integrar phenomenon over all distances and solid angles. We chose an iterative 
approach and used temperature as the convergence criterion. The iteration starts with the flame structure solver and 

an initial guess of = ~ea(T 4 -T*)/t to calculate velocities, gas and soot concentrations, density, and flame tempera- 
ture. The resulting temperature is used for calculating V-£ ra ^ from the YIX solver. The resulting V 'q Tac [ is substituted 
back into the structure solver to update (by averaging the two latest runs) velocities, concentrations, density, and 
temperature. This process is then repeated until a convergence on the temperature is accomplished between two 
subsequent iterations ( AT/ T< 0.002). The turbulence-radiation interactions, which have been shown to be significant 
by Gore et al. (ref. 27), will be included in the future. 

To evaluate spectral radiative properties, the absorption coefficient for soot aggregate is calculated from the Rayleigh 
solution (ref. 28) using the complex refractive index calculated from the Drude-Lorentz dispersion model and the 
parameter set of Habib and Vervisch (ref. 29). Scattering from soot particles is neglected at present to reduce computa- 
tional time, but will be included in the future. The exponential-wide-band model of Edwards (ref. 30) is used to 
calculate the absorption coefficient for C0 2 and H 2 0 gases. The most accurate way to evaluate the total (i.e., across the 
spectral range) fluxes and divergence of flux is first to calculate them for each subdivided spectral range and then 
calculate the quadratures, but this is inefficient due to the large number of spectral subdivisions necessary to account 
for the steep variation of gas absorption coefficients and the number of iterations between the flame structure and the 
radiation solver. To focus on the convergence of iterations, we will test the Planck ( a P ) and the Rossland ( a R ) mean 
absorption coefficients given as (ref. 31) 
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In the future, we will examine other more sophisticated approaches. 

Results and Discussions 

We first compare the radiation heat transfer results from the YIX method based on true spectral integration (1 to 20 
\im) and on both mean coefficients for a simulated turbulent ethylene diffusion flame. Fig. 1 shows a comparison of 
mixture spectral absorption coefficient against both means, and suggests that the Rossland mean should be reasona- 
bly accurate. Fig. 2 shows that the flux divergence based on the Rossland mean agrees quite well with that from 
integration. We thus base all subsequent "coupled" calculations on the Rossland mean coefficient. 

We next compare the soot formation models by Moss et al. and by Khan et al. in Fig. 3(a) for flame temperature and 
in Fig. 3(b) for soot volume fraction, against experimental data for a laminar flame (ref. 32). From these figures, it is 
unclear which model is superior. However, calculations based on the experimental temperature, which is higher 
than predicted values around the centerline, reveal that Khan’s model overpredicts the centerline volume fraction 
by a factor of 2, whereas Moss’s changes only slightly. We therefore favor Moss’s model at present. When we applied 
Moss's model with the same set of coefficients for the same fuel at different flow rates and compared with data, we 
found that the model is relatively insensitive to variations in fuel flow rate, and we will improve this in the future. 
Fig. 3(a) shows that the inclusion of radiative heat transfer provides a more accurate and correct trend in temperature 
profiles, when compared to those under adiabatic. Fig. 3(c) shows that the temperature converges smoothly in less 
than 10 iterations. 

We then tested the model on normal gravity (1-g) and microgravity (0-g) laminar acetylene flames, for which we 
have measured soot volume fraction data for comparison. Fig. 4 shows comparisons between 1-g and 0-g for (a) 
measured soot volume fraction, (b) predicted soot volume fraction, and (c) predicted temperature. The dotted curves 
in temperature map are predicted flow streaklines, which also indicate the paths traversed by soot particles. Results 
in (b) and (c) were calculated using soot formation modeling coefficients slightly different from those for ethylene in 
Fig. (3). It is encouraging that the soot formation model is capable of accurately predicting the overall level and shape 
of the volume fraction distribution, especially since this is a smoking flame. It is also encouraging that the same 1-g 
soot modeling coefficients can be used to predict 0-g volume fraction with reasonable accuracy. Although the model 
predicts a wider soot shell than that measured for both 1-g and 0-g, it is expected that further adjustments in soot 
modeling coefficients will yield even better predictions. 

Finally, we tested the model on turbulent ethylene flames, for which we should have measured data for comparison 
within a month or two. Fig. 5 shows a comparison of two predictions, one using the YIX solver and the other using 
the simple sink term model, against measured data for (a) volume fraction, and (b) temperature. Here, although it 
appears that the volume fraction results from the YIX solver are more accurate, the temperature results are roughly 
the same. Fig. 6 shows a comparison of 1-g and 0-g predictions for (a) volume fraction, and (b) temperature, using the 
same set of soot modeling coefficients and the YIX solver. Since this turbulent flame is momentum-dominant, the 
difference in flame shapes between 1-g and 0-g is not as substantial as that for the laminar flame in Fig. 4, which is 
more buoyancy-dominant. 

In our future work, we will develop more accurate means of evaluating total radiative quantities, such as the wide- 
band model used by Song and Viskanta (ref. 33). Their derivation of the "mean” radiative transfer equation to include 
turbulence-radiation effects is also worth noting. To investigate turbulence-radiation effects, we must also study the 
implementation of joint pdf (refs. 27 and 34) for correlating enthalpy, temperature, and mixture fraction. We will 
look for developments in turbulence-buoyancy interactions (ref. 15), which should be influential on 0-g flames. 
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Fig. 1 Absorption coefficients for a soot and gases 
mixture at 1400K, with 0.1 mole fraction for CO 2 and 
H 2 O, and a soot volume fraction of lOA 


Fig. 2 The divergence of flux contours for a simulated 
ethylene diffusion flame ( D = 0.58 mm, Q = 3.96 cm 3 / sec, 
and Re = 536). 
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(a) Flame Temperature 


(b) Soot Volume Fraction (c) Convergence of Temperature 


Fig. 3 Comparisons of model predictions to experimental data, and the trend of convergence of temperature during 
iterations between the structure and the radiation solver, for a laminar ethylene co-flow diffusion flame (ref. 32). 
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Fig. 4 Comparisons of model predictions to experimental data for a laminar acetylene jet diffusion flame. 


(a) Measured Volume Fraction 








Fig. 6 Predictions for 1-g and 0-g turbulent ethylene jet diffusion flames (D = 3 mm, V = 52 m/ s, and = 9.62x1 0 3 ). 
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